Charge and orbital order in half-doped manganites 
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An explanation is given for the charge order, orbital order and insulating state observed in half- 
doped manganese oxides, such as Ndi/2Sri/2Mn03. The competition between the kinetic energy of 
the electrons and the magnetic exchange energy drives the formation of effectively one-dimensional 
ferromagnetic zig-zag chains. Due to a topological phase factor in the hopping, the chains are intrin- 
sically insulating and orbital-ordered. Most surprisingly, the strong Coulomb interaction between 
electrons on the same Mn-ion leads to the experimentally observed charge ordering. For doping 
less than 1/2 the system is unstable towards phase separation into a ferromagnetic metallic and 
charge-ordered insulating phase. 



Manganese oxides with the general composition 
Ri_a:Aa;Mn03 (wheie R and A are rare- and alkaline- 
earth ions, respectively) have attracted considerable at- 
tention because of their unusual magnetic and elec- 
tronic properties. In some of these materials, metal- 
insulator transitions can be observed where both conduc- 
tivity and magnetization change markedly. The x = 
and X ~ \ end members of the Ri-j^AjjMnOs fam- 
ily are insulating and antiferromagnetic (AF) with the 
Mn-ion in the Mn'^+ and Mn'*+ state, respectively. For 
intermediate x, the average Mn- valence is non-integer 
and the material is generally metallic or semiconduct- 
ing. Most of the perovskite manganites show a ferro- 
magnetic (FM) groundstate when the holes are opti- 
mally doped (usually 0.2 < a; < 0.5) and anisotropic 
antiferromagnetic (AFM) phases for x > 0.5. The 
half-doped manganites, with x = ^, are very particu- 
lar. Magnetically these systems form FM zig-zag chains 
that are coupled AFM (see Fig.|l|) at low temperatures, 
the so-called magnetic CE-phase [|l|. The groundstate 
is, moreover, an orbital-ordered and charge-ordered in- 
sulator. This behavior is generic and is experimen- 
tally observed in Ndi/aSri/aMnOa ^j^, Pri/aSri/aMnOa 
i, Pri/2Cai/2Mn03 §, Lai/aCai/aMnOa §,0, 
Ndi/2Cai/2Mn03 [g[ and in the half-doped layered man- 
ganite Lai/2Sr3/2Mn04 ||]. The insulating charge- 
ordered state can be transformed into a metallic FM state 
by application of a external magnetic field, a transition 
that is accompanied by a change in resistivity of several 
orders of magnitude HJl^l • 

The occurrence of charge order, orbital order and large 
magneto-resistance in the half-doped manganites is ex- 
perimentally well established. Theoretically, however, 
the nature of the charge ordering at x = ^ and the origin 
of the unconventional zig-zag magnetic structure remain 
unclear. We address these points here and show that {i) 
the insulating CE-phase results from a particular order- 



ing of orbitals, (ii) this state is stable only in a narrow 
concentration range around the commensurate value and 
that {Hi) the zig-zag chains are intrinsically charge or- 
dered due to on-site Hubbard correlations. 

The competition between kinetic energy (double- 
exchange) and superexchange between the manganese 
core-spins, combined with the orbital degeneracy leads 
to the formation of the magnetic CE-phase, for the same 
reasons as it leads to the anisotropic magnetic phases in 
the highly doped manganites with x > ^ |T^. The CE- 
phase is found to be stable close to x = i; for x > 0.57 
the C-phase is stable and for a; < ^ we find phase sepa- 
ration between the CE-phase and FM-phase. 

As in the double-exchange framework electrons can 
only hop between sites with FM aligned core-spins, in 
the CE-phase only hopping processes within the zig-zag 
chains are possible, rendering the system one-dimensional 
for low-energy charge fluctuations. The essential obser- 
vation is that an electron that passes a corner-site of the 
zig-zag chain, acquires a phase that depends on the or- 
bital through which it passes. This leads to an effec- 
tive dimerization that splits the bands and opens a gap 
at the Fermi-surface. The gap is very robust as it is a 
consequence of the staggered phase-factor that is itself 
fully determined by the topology of the system. At the 
same time not all orbitals are fully occupied, leading to 
an orbital-polarized groundstate. Our most surprising 
observation is that the experimentally observed charge 
order is directly obtained from the degenerate double- 
exchange model when the Coulomb interaction (the Hub- 
bard U) between electrons in different orbitals, but on the 
same site is included. This can be understood from the 
fact that in the band-picture on the corner-sites both 
orthogonal orbitals are partially occupied, but on the 
bridge-site only one orbital is partially filled. The on- 
site Coulomb interaction acts therefore differently on the 
corner- and bridge-sites: charge is pushed away from the 
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FIG. 1. View of the CE-phase in the x-y plane. We choose 
our basis orbitals such that the gray lobes of the shown or- 
bitals have a negative sign. The dots at the bridge-sites rep- 
resent a charge-surplus. 



effectively correlated corner-sites to the effectively uncor- 
related bridge-sites. 

The theoretical explanation for the ferromagnetic 
metallic state in doped manganites was already devel- 
oped in 50's and 60's |12-1J]. The number of d electrons 
per manganese site is 4 — a; . Three electrons occupy the 
low lying t2g orbitals, having parallel spins. The other 
electrons occupy the eg orbitals with their spin parallel 
to the t2g spin because of the large ferromagnetic Hund's 
rule exchange Jh- In the double-exchange framework 
electrons can only hop from one Mn-site to a neighbor- 
ing one, when the t2g spins of these sites are aligned fer- 
romagnetically because otherwise, in the classical treat- 
ment of the t2g core-spin the electron would have to over- 
come an energy barrier proportional to Jh- 

On the other hand, the conventional superexchange fa- 
vors AF alignment the spins. This leads to an intricate 
competition between the kinetic energy and the superex- 
change [|lj| that is amplified by the two-fold degeneracy 
of the eg levels in a cubic crystal. As a basis for the Cg 
wavefunctions we take the x^ — y^ (|x)) &nd 3z^ — r^ (|C)) 
orbitals. The hopping of the electrons between neighbor- 
ing Mn-sites depends strongly on the kind of the orbitals 
involved and on the direction of the bond. The Hamilto- 
nian for the kinetic energy is 



Ht = iSo 



,/5,(u)r 



^4,/3 <aCj,/5, 



(1) 



where c^^^ creates an electron on site i in orbital a 
and (ij)r denotes a nearest neighbor pair along the F- 
dircction (F = x, y, z) and t is the hopping integral. The 
hopping is constrained to nearest neighbor Mn-ions that 
have FM aligned t2g core-spins (canted phases are not 
important in the doping range we consider here.) The 
matrices M^j^ and M^ a can be found by applying the 
cubic symmetry operations on M^^, where M|^> ,,> = 1 
and all other matrix elements are zero. This means phys- 
ically that along the z— direction, for instance, electrons 
can only hop between |^) orbitals and that all other hop- 
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FIG. 2. Top: topology of the interactions in a zig-zag chain, 
where ii = i/2, t2 — i-\/3/2, and U is the Coulomb interaction 
between electrons on the same site. Bottom left: electron 
dispersion in the zig-zag chain of the CE-phase for U — and 
electron dispersion in a straight chain, as in C-phase (dashed 
line). Bottom right: total energy per site for the CE-, C- 
and FM-phase for t/J = 5. The Maxwell-construction in the 
phase separated (PS) region is shown by the thin dashed line. 



pings (involving at least one |x) orbital) are zero by sym- 
metry, whereas in the x— and y— direction all hoppings 
are allowed (and related by symmetry). The strong spa- 
tial anisotropy of the hopping, combined with the com- 
peting kinetic and superexchange energy J can, depend- 
ing on the eg bandwidth and J, give rise to low dimen- 
sional spin-structures that, for instance, optimize kinetic 
energy by forming FM chains and optimize magnetic en- 
ergy by a AFM coupling of these chains (C-phase), or 
AFM coupled FM planes ( A-phase) (ll) . 

In the homogenous FM state at x = ^ the magnetic 
energy per site is E^^^ = 3J and the kinetic energy 
per site E^-^ = ~t, if the t2g spin is treated as a clas- 
sical spin. Here the magnetic energy is found by de- 
termining the number of AFM bonds, and the kinetic 
energy from the filling of the bands, as is described in 
Ref. |ll|. For AFM coupled FM chains, the C-phase, 
■^Sn + ^Siag = -0.6366i + J, SO that the C-phase is 
stable with respect to the FM-phase for J > 0.1816i. 
A third possibility is for the spins to form FM zig-za; 
chains that are coupled AFM, the CE-phase ||l|, see Fig 
The CE-phase is always higher in energy than the C- 
phase, except close to a: = i. For the half-doped system 
Egf + E^^g = -0.695i + J. It is easily shown that 
the A-phase is higher in energy. Thus for the half-doped 
manganites the experimentally observed CE-phase is in- 
deed lower in energy than the C-phase, and lower than 
the FM-phase when J > 0.1524t. 

Let us address our main findings that the CE-phase is 
charge-ordered, orbital-ordered and insulating, and come 
back to the stability of the CE-phase at the end of the 
paper. As shown in Fig.Q^the CE-phase contains two geo- 
metrically inequivalent sites, so called bridge- and corner- 
sites. Note that our specific choice of basis-orbitals as 
shown in Fig.n^ is motivated by the convenience of this ba- 
sis for the calculations. The expectation value of actual 



observables is, of course, independent from the choice 
of basis Wannier orbitals. The topology of the electron 
hopping integrals is shown in Fig.0. The crucial obser- 
vation is that an electron that hops from one bridge-site 
to another bridge-site via a |x) corner-orbital obtains a 
phase-factor —1, while if the hopping takes place via a 
10 corner-orbital, the phase-factor is -1-1. 

The bands are obtained by the solution of a 3x3 matrix. 
There are two bands with energy e± = ±ty/2 + cosk, 
where k is the wave vector, and two nondispersive bands 
at zero energy, see FigJ3. W^. At a; = ^ the e_ band 
is fully occupied, and all other bands are empty. The 
system is insulating as the occupied and empty bands 
are split by a gap A = t. In the C- and CE-phase the 
magnetic energies are equal, but the opening of the gap 
at the Fermi-energy in the CE-phase lowers its energy if 
the bands are half filled. This mechanism is equivalent 
to the situation in the lattice-Peierls problem, where the 
opening of a gap stabilizes a groundstate with a lattice 
deformation. 

In the insulating state the average occupancy of the 
bridge-site 3x^ — r^ {3y^ — r^) orbital is ^ while the or- 
thogonal bridge-orbital is empty. The average occupancy 
of the corner-site is also \, with the ratio between |x) 
and 10 occupancy of \/3 : 1. The system is therefore 
orbital-ordered, but charge is homogeneously distributed 
between corner and bridge-sites. 

This changes drastically when also the Coulomb inter- 
action U between electrons on the same site is taken into 
account. The Hamiltonian is 



Hij — UT^aKP.i ni^aTii^ii, 



(2) 



where rii^a = Cia'^i,a- For the tg electrons in the mangan- 
ites U K, lOf, so that the system is strongly correlated. 
The correlations, however, have a very different effect on 
the corner and bridge-sites. On the bridge-site, one or- 
bital is always empty, so that the Coulomb repulsion is 
ineffective, whereas on the corner-sites both orbitals are 
partially occupied. The consequence is that charge is 
pushed away from the correlated corner-sites to the ef- 
fectively uncorrelated bridge-sites, causing charge order. 
We treat the the full Hamiltonian H = Ht + Hjj for 
a; = i with three different methods: exact diagonaliza- 
tion (ED), in mean- field (MF) and with the Gutzwiller- 
projection (GP) ||l3]. In the ED calculation we consider 
a ring of 12 and 14 sites (or 18 and 21 orbitals, respec- 
tively) with the topology as in FigJ| . In Fig.y the charge- 
disproportionation 5 as a function of U/t is shown, where 
6 is defined as {uB — nA}, where {ub) {{nA)) is the expec- 
tation value for finding an electron on a bridge- (corner- 
) site. The results for the 12 and 14-site cluster differ 
by less than 2% and finite size effects are therefore very 
small. For U/t = also (5 = 0, as explained above, and 
for finite Coulomb interaction 5 monotonously increases 
with U . For U/t = 10, (5 = 0.171 and 5 reaches its maxi- 
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FIG. 3. Charge-disproportionation as a function of 
Coulomb interaction and doping (inset). The full line is ob- 
tained from exact diagonalization (ED) of a 14 site cluster. 
The dashed lines for small and large U are obtained by the 
mean-field and Gutzwiller approximations, respectively. The 
dots in the inset are the ED results for a 12 and 14 site cluster 
for U -^ CO. 

mum value for U -^ oo. In a ED calculation for a 16-site 
cluster with U -^ oo, we find that 5^ = 0.185. 

For small values of U/t a MF treatment is expected 
to be very reliable because the charge fluctuations are 
gapped. We decouple the quartic term in Eq.2 and deter- 
mine the charge densities in the orbitals self-consistently. 
The results for small U/t are shown with a dashed line 
in FigJS and they agree very well with the ED results for 
U < Is.. In the GP, valid for C/ ^ oo, we introduce con- 
strained electrons on the correlated corner-site. For the 
Ix) orbital we introduce Xi = {1 — nl)xi and for the \C/) 
orbital Zi — {1 — n^)zi, where Xi (zi) denotes an electron 
operator on corner-site i, with corresponding density nf 
{nl) and Xi (zi) a constrained electron operator. The 
Hamiltonian is 



Hgp = SieB 4 [ti{zi+i + Zi^i) 

+ t2{Xi+i - Xi_i) ] + h.c, 



(3) 



which physically means that an electron can only hop to 
a |x)-orbital on a corner-site when there is no electron 
in the |C)-orbital on that site and vice versa, so that a 
corner-site can never be occupied by more than one elec- 
tron, as is required in the [/-infinity limit. We decouple 
the quartic terms in Hgp as chj^Ziuf 



W+l^M 



{Cij^iZi)n^ . After a Fourier transform we find 

+ iik{c\zk + zlck) + hkiclxk + xlck), 



(4) 



with ex = -4iiEfcC0sfc (c[zfc), ez = 4it2SfeSinfc (cj^Xfc), 
tik — 2ti{l — n^)cosk and t2k — —2it2{l — n^)smk. 
We obtain the solution of this system of equations by 
iteration, and find the charge-disproportionation S^^ = 
0.191, which is represented by the dashed line in Fig.H 
for large U/t, and agrees well with the ED results. 

The on-site Coulomb interaction causes a charge- 
surplus on the bridge-sites in the x-y plane. In the 



CE-structure the zig-zag chains are stacked AFM along 
the z-direction, which imphes that above each bridge- 
site there is another bridge-site in the next plane. So 
we find that the charges actually accumulate on sheets 
formed by the bridge-sites along the z-direction. This 
is in remarkable agreement with experiment and at the 
same time excludes the possibility that the charge-order 
is driven by longer range Coulomb interactions because 
the Madelung-sum is always minimized for a rocksalt- 
type charge-order. Similar physics may also apply to the 
situation at x > 0.5, e.g. for the stripe phases observed 
in Ref. §. 

Fig.y shows that the charge disproportionation is 
strongly doping dependent. For x > ^ the holes that are 
doped into the lower e_-band efficiently suppress charge- 
order. In this doping range, however, the CE-phase be- 
comes unstable with respect to the C-phase. In Fig.|| is 
shown that the kinetic energy of the C-phase is lower for 
X > 0.57. For x < -^ the energy per site of the CE-phase, 
EcE, is constant because the extra electrons are doped in 
the non-dispersive bands at zero energy, which causes a 
kink oi EcE at a; = i. For lower hole-doping (higher elec- 
tron concentration) the homogeneous FM-phase is lower 
in energy, as is expected. The stability of the phases can 
be tested by performing a Maxwell-construction. Due to 
the kink in Ece at x = ^j one finds that for x < ^ 
the FM-phase and CE-phase coexist. The doping region 
where this phase separation occurs depends strongly on 
the actual value of the electron bandwidth and magnetic 
exchange. We conclude that there is a strong asymmetry 
for doping lower and higher than i, as phase separation 
between the CE- and FM-phases is only present on the 
lower doping side. Inclusion of the correlations beyond 
the mean-field approximation and longer range Coulomb 
interactions, which we did not consider in the calcula- 
tion of total energies, may modify these results. On the 
basis of the ED results we believe, however, that U has 
only a small effect on the total energies. A more rele- 
vant contribution to the total energy comes from lattice 
deformations, as it is obvious that in the Jahn- Teller dis- 
torted C- and, to a lesser extend, CE-phases [Q there is 
a large gain in lattice energy with respect to the undis- 
torted FM-phase. It can be shown that the lattice contri- 
butions to the groundstate energy can, to a large extend, 
be described by using effective values of J [1^], so that 
our general conclusions about phase separation for x < -^ 
and about the asymmetry of the phase diagram for x < ■^ 
and X > ■^ remain valid. Recently this phase separation 
into a FM metallic and a charge ordered AFM insulating 
phase was observed experimentally p9[ | . 

In summary, we have given an explanation for the 
charge order, orbital order and insulating state of the 
half-doped manganese oxides. The CE-phase with one- 
dimensional ferromagnetic zig-zag chains is stable for 
X = -^ and due to staggered phase factors in the hopping, 
the chains are insulating and orbital-ordered. The strik- 



ing feature of our model is that the strong Coulomb inter- 
action between electrons on the same Mn-site leads to the 
experimentally observed charge ordering. In a magnetic 
field we expect the chains to be unstable with respect 
to the ferromagnetic metallic state. This might offer a 
likely explanation for the large magneto-resistance at the 
metal-insulator transition for the half-doped manganites. 
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